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We analyze two-component spatial optical vortex solitons supported by parametric wave mixing 
q;^", processes in a nonlinear bulk medium. We study two distinct cases of such localised waves, namely, 

' parametric vortex solitons due to phase-matched second-harmonic generation in an optical medium 



(N 



with competing quadratic and cubic nonlinear response, and vortex solitons in the presence of third- 
harmonic generation in a cubic medium. We find, analytically and numerically, the structure of 
two-component vortex solitons, and also investigate modulational instability of their plane-wave 



' background. In particular, we predict and analyze in detail novel types of vortex solitons, a 'halo- 

■ vortex', consisting of a two-component vortex core surrounded by a bright ring of its harmonic field, 

o . 

, and a 'ring-vortex' soliton which is a vortex in a harmonic field that guides a ring- like localized 

' mode of the fundamental-frequency field. 

g\ 

^ . I. INTRODUCTION 

» 
> 

. . . 

O |. An optical vortex soliton appears as a stationary self-trapped beam in a self-defocusing optical medium that carries a 
IL^ phase singularity on an electromagnetic field, so that the beam intensity vanishes at a certain point, and the field phase 
changes by 27rm (m being integer) along any closed loop around the zero-intensity point. If such an object is created 
> in a linear bulk medium |l|,|^, it preserves the singularity but expands due to diffraction. However, in a nonlinear 
?H medium, the diffraction-induced expansion of the vortex core can be compensated for by a nonlinearity- induced change 
. . . in the refractive index of a nonlinear medium, thereby creating a stationary self-trapped structure, an optical vortex 
soliton. Such nonlinear localized waves carrying a singularity were first introduced as stationary solutions of the 
nonlinear Schrodinger (NLS) equation in the pioneering paper by Ginzburg and Pitaevsky 1^ to describe topological 
excitations in superfiuids, but the same objects appear in many other fields including nonlinear optics 

The parametric interactions may provide an efficient way of vortex transformation. In particular, by mixing waves 
of different frequencies, one can change the vortex topological charge m and even the vortex polarization. Recently, the 
first experimental results on the vortex generation in the presence of two- wave parametric mixing have been reported in 
nonlinear optics, including the second- harmonic generation (SHG) and more general types of frequency conversion 
1^ and sum-frequency mixing where the generation of higher-order {\m\ > 1) linear vortices in the case of negligible 
spatial walk-off between harmonics was demonstrated. 

To the best of our knowledge, no theory of parametric optical vortices in the presence of both diffraction and 
nonlinearity has been developed so far. In a nonlinear regime, an interplay between diffraction and parametric 
coupling of the harmonic fields is expected to lead to the formation of stationary structures - parametric vortex solitons 
- supported by three- or four-wave mixing between the phase-matched waves of different frequencies. Stability of such 
multi- frequency vortex solitons is a key issue. For example, in the problem of SHG in a diffractive bulk medium, 
vortex solitons are expected to be unstable due to parametric modulational instability of the two-wave background 
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field |10|. Recently, it has been suggested that taking into account a weak defocusing cubic nonlinearity one 
can eliminate the development of parametric modulational instability allowing stable dark solitons to exist. Some 
examples of stable two- wave parametric dark solitons have been presented in Ref. |ll|] , and it has been pointed out 
that, in the problem of SHG, a stable vortex soliton of the lowest possible charge (|m| — 1) can exist describing a 
27r-phase twist of the fundamental wave and 47r-phase twist in the second-harmonic field. 

In the present paper we suggest a general approach to the analysis of multi- component vortex solitons resulting from 
parametric wave mixing. The general theory is then developed in detail in the no-walkoff case for two examples: (i) 
parametric interaction of the first and second harmonics in a medium with competing quadratic and cubic nonlinearity, 
and (ii) parametric interaction between the first and third harmonics in a medium with a cubic nonlinear response. 
In both the cases we find different classes of vortex solitons as (2-1-1)- dimensional dark solitons of circular symmetry 
carrying a phase singularity, and investigate their stability to propagation and modulational stability of the supporting 
two-component background waves. 

The paper is organized as follows. In Sec. II we briefly present two models of parametric wave interaction that 
describe a phase-matched coupling between the fundamental frequency mode and its harmonic field, in the case of 
phase-matched wave mixing and no walk-off. The further analysis of the asymptotic structure of stationary localized 
solutions for parametric vortex solitons is rather general, and it is presented in Sec. Ill for both the models. Section 
IV is devoted to the analysis of vortex solitons in the model of competing nonlincarities. We find numerically the 
profiles of two-component vortex solitons and investigate their stability to propagation. In particular, we reveal the 
existence of novel classes of dark-soliton solutions of radial symmetry, including a ring-vortex soliton, that consists 
of a vortex core in the harmonic field surrounded by a bright ring of its fundamental frequency, and a halo-vortex, a 
two-wave vortex soliton with nonmonotonic tails. The corresponding results are also obtained for the problem of the 
third-harmonic generation in Sec. V. Finally, Sec. VI gives the summary of our results and briefly discusses some 
related issues including the comments on experimental verifications and a link with other problems. 



II. MODELS OF TWO- WAVE PARAMETRIC INTERACTION 



A. Competing Nonlinearities 



First, we consider the model of competing quadratic and cubic nonlinearities introduced earlier for the (1+1)- 
dimensional case in Ref. (l2j and recently generalized to the case of (2+l)-dimensional bright solitons of radial 
symmetry in a bulk medium p3[ . We assume that a beam of a fundamental harmonic (FH) with the frequency 
is launched into a medium possessing combined quadratic [or x^^"*] ^^^id cubic [or x''^^] nonlinear response under the 
condition of phase-matched SHG. The FH beam generates a second harmonic (SH) wave, and such a two-wave mixing 
process in a bulk medium is described by a system of two coupled nonlinear equations. 
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(\E2\^ + p\Ei\^)e2^0, 



where Ei and E2 are the complex amplitude envelopes of FH (wi = uj) and SH {lj2 — 2a;) waves, respectively; 
ki — k{u!) and k2 = k{2uj) are the corresponding wave numbers; Afc = (2fci — ^2) is the wave-vector mismatch 
between the harmonics, p (which we take p = 2) is the cross-phase- modulation coefficient, and the coefficients x'^' 
and x^^^ are proportional to the second- and third-order susceptibility tensor elements and they characterize the 
combined nonlinear response of an optical medium. 

Adopting a similar set of scaling transformations as in Ref. [ p^ , we measure the transverse coordinates in the units 
of the beam radius Rq, and the propagation coordinate, in the units of the beam diffraction length = 2fcii?Q. 
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Then, applying the transformations 

where the parameter /3 stands for the nonhnearity-induced change of the beam propagation constant and A — 
2fcii?^Afc, we obtain a system of normalized equations for u and 



(2) 



du 2 * /I^P i2\ 

i-g^ + sV _i^u — u + wu '^Xy-^+P\w\ J" = 0, 

• ^2 /r, I l2 I i2n 

la— h sV^it; — au) + — + x(2(t|w| + p|u| )w = 0, 

where a = (2/3 + A)cr//3, s = sign/?, and the coordinates are rescaled as follows z z/j3 and {x,y) {x,y)/y 
For the spatial beam propagation we take a — 2. Parameter x describes a competition between quadratic and cubic 
nonlinearities, and it is defined as 



^ '^167rL^2^2 [^(2)]2- W 



Stationary solutions are then described by Eqs. with the z-derivatives omitted. To look for radially symmetric 
solutions carrying a phase singularity, we use the polar coordinates r = \J x^ +2/^, — tan^^(a;/y), and make the 
following substitutions, 

M(r, 4>) = [/(r)e''"'^, u;(r, 0) = W^(r)e^^'"'^, (4) 

where U (r) and W(r) are real functions and, for parametric interaction between the fundamental and second har- 
monics, iV = 2 whereas m is an integer number that characterises the vortex charge. 

Substituting Eqs. (^) into Eqs. (|2|), we obtain 

d^U IdU m^U^ -0 
dr"^ r dr dU ' 

d^W 1 dW m^N'^W^ -0 
dr^ r dr dW ' 

where the function F has the meaning of an effective potential, and it is defined as 

F = Fi(U, W) = --U^ + -U'^W - -W^ 
2 2 2 



B. Third-Harmonic Generation 



A similar type of two-wave parametric interaction occurs under the condition of the third-harmonic generation 
(THG). Bright and dark solitary waves in a waveguide geometry (i.e. with one transverse dimension) have been 
analyzed in Ref. p4| . In this case, the parametric interaction occurs between the fundamental beam (uji — uj) and 
its third harmonic(a;3 = 3a;), and the corresponding physical model of the parametric wave mixing in a bulk can be 
described by a system of two coupled equations. 
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where Ei and E3 are the slowly varying envelopes of the first and third harmonic fields, respectively, with corresponding 
wave numbers ki = k{uj) and = fc(3a;); A/c = 3fci — k^ is the wave- vector mismatch between the harmonics and 
X — {iTTUj"^ / c^)\x^'^^\ is the nonlinearity parameter, which is assumed here to be always positive, whereas x^^'' < 0. 

We follow a normalisation procedure similar to that used above for the competing nonlinearity model. Again, the 
transverse coordinate is measured in units of the beam width i?o and the propagation coordinate, in units of the 
diffraction length Rd — 2fcii?Q. Using the transformations of Rcf. iQ 



El = i y/f3/3^kiRlx] e'^'u{x, y, z), 
E2 = [V^I\lkiRlxj e^^^P+^>M^,y,z), 



the physical equations can be written in the following normalised form [cf. Eqs 

.du 



+ sV^u - u - -u*^w - s ( ^ + 2\wf I u = 0, 
oz 3 

dw 



(9) 



icr^ + sV^w -aw - -u^ - s(9|wP + 2\u\^)w 



where u and w are the normalised amplitudes of the fundamental harmonic field and its third harmonic, a — a(3P + 
A)//3, A = 2fcii?§Afc, s = sign/3, the transverse and propagation coordinates have been rescaled in terms of the 
nonlinearity-induced change of the propagation constant /3, 2 — > z/(3 and {x,y) {x , y) / ^/\(3\ , and, for spatial 
solitons, we take cr = 3. Importantly, everywhere below we consider only defocusing cubic nonlinearity searching for 
vortex-type solitary waves on a modulationally stable nonvanishing background. 

Stationary radially symmetric localized solutions of Eqs. (^) have the form (Q) with iV = 3, and they satisfy Eqs. 
with the potential F, this time defined as 

F ^ F2iU,W) ^ -^U^ - ^W^ 

^ ^ (10) 

-s(lu^W + + ^W^ + w^uA . 
\9 36 4 / 

Thus, in both the cases, stationary vortex-like structures are described by the same system of equations (|^) with 
different types of the potential F. This observation allows us to perform further analytical calculations in a rather 
general form, and, therefore, most of them are universal and can be applied to other models. 



III. GENERAL THEORY OF PARAMETRIC VORTEX SOLITONS 



A. Stationary Solutions 



Stationary radially symmetric solutions of Eqs. (||) [Eqs. (^)] are given by Eqs. ^ with the potential function F 
defined in Eq. (^) [Eq. (p^)] and N — 2 [iV = 3]. It is important to note that the parametric coupling between the 
modes brings several new features in the vortex structure and properties. Indeed, as follows from Eqs. (^) and (j^), a 
vortex with the charge m in the fundamental mode is always coupled to a vortex of the charge Nm {N = 2, 3) in the 
harmonic component. This makes parametric vortices very different from all types of vortex solitons analyzed earlier 
in the systems of two incoherently coupled NLS equations (see, e.g., Ref. and references therein). 
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B. Analysis of Vortex Asymptotics 



We are interested in the localized solutions supported by a two-component finite-amplitude background wave. For 
r — > 00, the background amplitudes (C/qj Wq) satisfy the coupled algebraic equations: 

dF dF 

which may have one or more nontrivial solutions. Importantly, due to the self-action effect we always have a special 
solution of the form (0, Wq), that corresponds to an excited harmonic field only. 

A vortex soliton is a localized nonlinear mode that asymptotically approaches the background (C/q, Wq) for r — > 00, 
but its intensity vanishes for r ^ to keep the terms ~ {ni^/r^)!! and ~ {m'^N'^ /r'^)W in Eqs. finite. This 
implies that we can find the vortex asymptotics in a rather general form. For r ^ 0, we look for solutions of Eqs. (|^) 
in the form: 



W = Wo 



2 



A A 



r2 



where (C/q, Wq) is a solution of Eqs. for the background amplitudes. Keeping in Eqs. (g) only the asymptotic 
terms up to the order of l/r^, we obtain, 



9 rd'^F\ , ( d^F \ 

where the index '0' stands for the values calculated &X\J — Vq and W — Wq. Solutions of the linear equations ( p^ ) 
for A and B can be easily found analytically; they define the asymptotics of the vortex solitons for different values of 
the vortex charge m in terms of the background amplitudes Uq and Wq defined by Eqs. (pj]). 

The analysis of the asymptotics gives us important information about the vortex structure. If both the products 
AUq and BWq are positive [see Eqs. (p^j, the vortex has a standard profile with the intensity in the core growing 
monotonically and always lower than the background intensity. However, if one of these products is negative, some- 
where across the vortex the intensity becomes higher than the asymptotic background intensity. That implies that 
the vortex core is surrounded by a bright ring of higher intensity. We call such structures 'halo-vortices '. In both the 
cases mentioned above, such vortex solitons may exist on a modulationally stable background, and some examples 
are given below in Sections IV and M. 



C. Vortex Soliton as a Waveguide 



The concept of light guiding light (see e.g., Ref. and references therein) is based on a simple observation that 
a spatial optical soliton (e.g., vortex) creates an effective optical waveguide in a nonlinear medium that can guide a 
wave of different frequency or polarization. It is clear that a vortex soliton creates a waveguide of radial symmetry 
which can guide a fundamental mode (no nodes) of the other wave. For the case of two incoherently coupled NLS 
equations describing two orthogonal polarizations, the guiding properties of vortex solitons have been analyzed by 
Haelterman and Sheppard | p5[ . The first demonstration of an optically written waveguide based on an optical vortex 
has been recently reported by Truscott et al However, the theory developed in Ref. (l^ is not valid for the case 
of the resonant interactions and parametrically coupled waves. Indeed, the parametric interaction forces the harmonic 
field to vanish for r ^ 0, trapping a singularity of the order of Nm. Therefore, a parametric vortex cannot guide a 
fundamental mode. To analyse the guiding properties of parametric vortex solitons, we note that Eqs. (^) with the 
potential F defined by Eqs. (||) and ( |l0[ ) have the solution {Uo — 0, Wq ^ 0). Therefore, we consider a vortex soliton 
created by a harmonic field W , with a stationary profile described by the nonlinear equation, 

-— + -— - + sa] W - jW"" = 0, 

dr'^ r dr \ r'^ I 
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where 7 = — 4sx, for the model (||), and 7 = 9, for the model (^. This equation always has a solution in the form 
of a vortex soliton with the charge Nm provided 7 > and sa < 0. Now, an eigenvalue equation for a linear mode 
guided by the vortex W{r) follows from the first equation of the system (||). Assuming U <€. max(W), we obtain, 



where 



1 dU 
r dr 



G{r) 



/d^F\ 
\dlp) 



sG(r) 



[/ = 0, 



(14) 



u=o, 



Equation ( p^ ) is a standard eigenvalue problem of the linear waveguide theory, and it can be studied analytically, 
e.g. by means of variational methods (see, e.g., Ref. |l8| and references therein). To make some analytical estimates, 
we present G(r) in an approximate form and obtain 



dr^ 



1 dU 
r dr 



-U -EU 



(L>2 



-U = 0, 



where E, C, and D are, in general, functions of a and 7. The parameters are chosen to provide the best approximation 
of the effective potential G{r). Using the standard variational method (or Ritz optimisation approach) and looking for 
a bifurcation of a linear mode taken in a trial form, /(r) = rexp (— Kr), we obtain an implicit expression to determine 
the mode cutoff a, 

E = (2(7^ - 1 - 2m^f /{mC^D^), 
which we analyse below for some particular cases. 



D. Modulational Instability 



Stability of the stationary vortex solitons described by the system (^) is an important issue. In general, the stability 
analysis of vortices in nonlinear models is a complicated and, generally speaking, unsolved problem. Instability can 
develop due to the presence of unstable eigenmodes localized near the vortex core and, in the one-dimensional case, 
this type of instability of dark solitons leads to the soliton motion, i.e. it is a drift instability (see, e.g., Ref. and 
references therein). Since moving vortices with nonzero minimum intensity (similar to grey solitons) do not exist, 
similar drift instability is not observed for vortices. The main instability which is usually associated with a vortex 
soliton originates from the instability of the nonlocalized background field. 

The analysis of modulational instability of the background field can be carried out in a general form. First, we 
write Eqs. (||) and Eqs. (^) in the form 



with T defined as 



du „9 dJ- 
i— + sV^u+— = 0, 
dz ou* 

la— + sV^w + - — = 0, 
dz aw* 



T ^Ti = + ^(u^w* + u*2w) - a\w\^ + 

x[\\ut + 2\w\^ + p\w\M''), 



(15) 



(16) 



for the model of competing nonlinearities, or 



T ^T2 = + ^(u^w* + u*^w) - a\w\^ + 



(17) 
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for the model of the third-harmonic generation. We look for stability of the background wave solution {Uq, Wq) defined 
by Eqs. (|ll|), and linearize Eqs. around this stationary solution substituting: 



u = Uo + ae"'-''+'"^' + he 

As a result, we obtain a system of linear equations for a, 6*, c, and d* leading to the characteristic equation: 



(18) 



4.*— O 4. 4.. 4. 

A * ^ 4* 4.._0 4* 

-^w.u* -^w.u -^w.w* -^w.w ^ 



= 0. 



Here 



{d"^ T / dn d'm)\(^u=Uo,w=Wo)i where m,n= {u,u* ,w,w*) {n ^ m), and for the a = f3 case 



\{u=Uo.w=Wo) 



Solving the characteristic equation with respect to O, we conclude the modulational instability analysis: purely real 
fl solutions for all positive |kp (with other parameters fixed) indicate a modulationally stable background for this 
fixed set of the parameters. Below we present the results of the modulational instability analysis for two cases of the 
parametric two-wave interaction. 



IV. COMPETING NONLINEARITIES 



Analysis of modulational instability for the system (^) has been briefly presented in Ref. ||Tl| . Below we repeat the 
main steps of that analysis for the completeness of this paper. Solutions of Eqs. (||) for background waves can be 
found by solving the coupled algebraic equations: 

12x< + 12< + Wo{a - 8 + 2x-') = 2x'\ 

ul = 4(1 - Wo)x-^ - ml 



for real J7o and Wq. There exist up to three such solutions with both amplitudes [/o and Wo being nonzero. Performing 
the analysis of modulational instability for each of the three solutions at s = ±1, we find that there exists only one 
modulationally stable mode. The parameter domains where such a solution exists are presented in Figs. [|(a,b). Note, 
that sign(sx) = signx'^'^^ and thus modulationally stable solutions exist only for x^'^^ < 0. Importantly, the amplitude 
of the modulationally stable background diverge in the limit x — > so that the stable background solution exists 
exclusively due to mutual action of quadratic and cubic nonlinearities. Other nonlinear modes are modulationally 
unstable in the whole domain of their existence and they are not presented in Figs. 0(a,b). Modulational stability 
in the limit of large negative x^^'' [e.g., for s — —1 and x > 0; see Fig. |^(a)], is not surprising because stable dark 
solitons are known to exist in a defocusing Kerr medium without quadratic nonlinearity. Here, we are interested in 
the case when the effective nonhnearity is predominantly quadratic, i.e. IxC^ol ~ Ix^^ol ^ 1- We found that this 
condition can only be satisfied for s = +1 where modulationally stable background waves of moderate amplitudes 
exist for relatively small values of negative x [see Fig. 1(b)]. 

Using the numerical relaxation technique, we have found that a continuous family of two-component vortex solitons 
exists in the whole region of the existence of modulationally stable background waves shown in Figs. |l|(a, b). Figures 
§(a) and |^(b) present an example of such a vortex soliton. 

Analysis of the vortex asymptotics demonstrates that halo- vortices can exist in Eqs. (||) only if s = — 1 and x > Oi 
in the narrow domain shown in Fig. |^. Both terms with B and B2 factors in the asymptotic expansion ( |T^ ) contribute 
to the formation of the halo (i.e. both BWq and B2W0 products are negative). 

Ring-vortex solitons can also exist in the model (||), see Figs. 4(a,b) and Figs. 5(a,b). Variational analysis allows 
us to find an approximate expression for the bifurcation curve where such solutions appear, [2 ~ a) — ^ ct/x- As 
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a increases, the maximum of the bright-ring amphtude approaches the value /7o of two-wave modulationally stable 
parametric plane waves. At the values of r where IJ approaches J7o, the second component, W ^ also approaches the 
corresponding plane wave amplitude Wq [see Fig. ib)]. 

Such ring-vortex solitons can be unstable due to modulational instability of the background wave Uq — 0, Wq — 
aj^X- For example, for s = —1, a > 0, modulational stability is defined by the condition (a — 2) > yja/x- The 
regions of the existence of modulationally stable one-component plane waves and ring-vortex solitons of Eqs. are 
presented in Fig. ||. 

Existence of two-component stable vortex solitons composed of parametrically coupled fields suggests that such 
vortices can be excited in the process of the harmonic generation. In Fig. 7 we present the numerical simulation results 
supporting this idea. We launch a mode of the fundamental frequency without a seeded second harmonic assuming the 
condition of phase-matching. The vortex soliton dynamics is simulated using a split-step beam propagation method 
(BPM). To solve the problem of the vortex phase geometry, we simulate a system of four vortices, arranged such that 
horizontally and vertically adjacent vortices are opposite in charge. Lines of equal phase are chosen to correspond to 
the lines of the force of an equivalent system of electrostatic point charges. With the periodic boundary conditions 
imposed by BPM, this configuration means that, in fact, an infinite array of vortices is simulated. Figure 7 shows the 
vortex generation by an input fundamental mode with single-charged vortices. Due to phase matching with the second 
harmonic, we observe a generation of double-charge vortices in the harmonic field (see the plots at z = 1) and then 
periodic oscillations of the two-component background and the vortex profiles near a stationary state corresponding 
to a lattice of two-component vortex solitons (see the plots at z = 10 as an example of such dynamics). 



V. THIRD-HARMONIC GENERATION 



Vortex solitons of Eqs. (|^) have fewer parameters in comparison with the parametric vortices described by Eqs. 
, and thus they can be analysed much more easily numerically. These vortex solitons are found in the whole region 
of the existence of modulationally stable plane waves, i.e. for a < ath ~ 14.509. Examples of such vortex solitons are 
shown in Figs. ||(a, b). 

A third-harmonic component of the vortex solitons has a nonmonotonic tail for 10.85 < a < 14.509, however it can 
only be called a halo vortex for the interval 11.26 < a < 14.509, where the absolute value of a local extremum in the 
structure of the vortex tail is greater than the corresponding plane wave background value Wq. The halo is becoming 
more pronounced as a ^ 14.509. An example of a halo-vortex soliton of Eqs. ^ is shown in Fig. ^. 

Ring-vortex solitons have also been found for the model (^), see Fig. 10. In this case, a variational analysis allows 
us to find an approximate analytical result for the bifurcation point (s = —1) where a ring-like mode is guided by the 
vortex: 

(jB 

abif = o— • 

For 7 = 9, (5 = 2, A^ = 3, and m = 1 this gives abif ~ 4.52, which agrees well with numerical data. We also find that, 
in general, coupled ring-vortex solitons exist for a > 4.5, and for each such value of a (except a = abif) there exist two 
different types of ring-vortex solitons (see Fig. p^ ). As the parameter a decreases, the maximum of the bright ring 
in the fundamental mode approaches the value of Uq of the two-wave modulationally stable background field. Again, 
as it has been observed for the model of competing nonlinearities, at values of r where U approaches Uq, the vortex 
component W deforms significantly approaching the corresponding plane- wave amplitude Wq. We note that all these 
ring-vortex solitons are modulationally stable, because, in the framework of Eqs. (|^), modulational instability does 
not occur for one-component plane wave solutions. 



VI. CONCLUDING REMARKS 



We have analyzed two-component vortex solitons supported by parametric wave mixing in a nonlinear optical 
medium. We have considered two classes of such vortex solitons. In the first case, we have studied the existence, 
structure, and stability of vortex solitons supported by phase-matched interaction between the fundamental and 
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second-harmonic waves in a quadratic medium, and the effect of the next-order cubic nonhnearity has been taken into 
account for suppressing modulational instabiUty of the supporting plane-wave background. In the second case, we 
have considered how the vortex parameters, structure, and stabihty are modified due to the process of third-harmonic 
generation when the phase-matched wave interaction generates a corresponding multi-charge vortex component in 
a harmonic field. In particular, we have predicted the so-called 'halo-vortex' consisting of a two-wave vortex core 
surrounded by a bright ring on a non- vanishing background. Additionally, we have analyzed the waveguiding properties 
of a vortex soliton in the case when it guides a harmonic field due to a phase-matched parametric interaction. A 
rigorous analysis of the stability of these parametric vortex solitons is still an open problem, as well as the effect of 
walk-off on the vortex existence and stability. 

As for experimental verifications of the vortex solitons described above, we would like to mention that, at least in 
the low-intensity regime, parametric vortices have already been observed in nonlinear optics. A possibility of SHG by 
a beam with a vortex was first mentioned and experimentally verified in Ref. Q, where a vortex of the topological 
charge m = 2 was found in the second-harmonic wave when the fundamental wave contained a vortex of the topological 
charge m — 1. The similar results on SHG have been presented by Dholakia et al in Ref. , whereas more complicated 
processes of sum- frequency mixing with beams carrying phase singularities were reported by Berzanskis et al . It is 
worth noticing that in all of those observations the different harmonics experienced noticeable walk-off that makes the 
stationary structures difhcult to observe, also introducing novel features in the vortex dynamics. In particular, for a 
collinear type I phase-matched SHG with an input beam carrying a single-charge vortex, Matijosius et al observed 
two intensity zeroes in the second-harmonic field with the separation of two SHG vortices due to walk-off. Thus, we 
can expect that stationary two-component vortex solitons discussed above can be observed in typical upconversion 
experiments when a high-intensity beam undergoes frequency doubling simultaneously with the creation of a phase 
singularity produced by a phase mask at the input, similar to the experiments mentioned above which were performed 
at moderate powers. Stability of those vortex solitons requires small (or zero) walk-off and a small defocusing Kerr 
nonlinearity of an optical material at both (or at least fundamental wave) frequencies. 

Additionally, we would like to mention that the parametrically coupled equations of competing nonlinearities, similar 



to Eqs. (^ analysed above, have been recently introduced by Heinzcn et al |20| to describe the dynamics of coupled 
atomic and molecular Bose-Einstein condensates, leading to a kind of "super-chemistry" in which the formation of 
molecules is a controlled parametric quantum process. In spite of the fact that both atomic and molecular condensates 
should be considered in a trapping external potential |21| , many of the features of the coupled stationary states, 
including all the types of the vortex states introduced above, are expected to exist in the model of atom-molecular 
condensates as well, providing a much broader view of the phenomenology of parametric vortex solitons. 

At last, we expect that the concept of the two-component parametric vortices, generated and supported by the third- 
harmonic generation process, can be important in the so-called third-harmonic microscopy (see, e.g., [Q) where an 
image is rendered using a series of cross-sectional images produced by third-harmonic generation within the specimen. 
Vortices can then be formed due to the development of caustics Q| in the reflected harmonic field, indicating the 
regions of highly concentrated inhomogeneities. This technique is based on the fact that the nonlinear susceptibility 
of solid media vary over many orders of magnitude, compared with linear refractive index changes that vary in a 
relatively small range. 
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FIG. 1. Existence domains for the modulationally stable background modes of the system 
s — +1 and X < 0- 
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FIG. 2. An example of a two-component vortex soliton supported by competing nonlinearity (a 



-2.5, X 



-0.1, and 



s = -1-1). In (a) only the vortex of the fundamental harmonic field is shown. 
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FIG. 3. Region of the existence of halo- vortices in the model of competing nonlinearities, s = —1. 



0.8 



O.Q 



♦ 
♦ 

* 
1 
' 

/ \* 
/ ' \ 




/ ' \ 
f ' \ 
1 ' \ 
1 1 \ 
1 • \ 
la % 
f ' \ 
f * \ 

1 * \ 

/ \u 


(a) 








FIG. 4. Examples of a ring- vortex soliton (s = — 1, x = 1) for (a) a = 1.4 and (b) q = 3. Solid - FH wave, dashed - SH wave. 
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FIG. 6. Region of existence for ring-vortex solitons at s = — 1 and x > 0. Left curve, (2 — a) = \J ajx-, is defined by a 
bifurcation where a bright component appears. Middle curve is a critical threshold that divides modulationally stable and 
unstable solutions. Right-hand side curve is the boundary for the existence of two-component parametric plane waves. 
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FIG. 7. Generation of two-component parametric vortex solitons by the fundamental vortices. Profiles of the fundamental 
(upper row) and second-harmonic (lower row) fields are shown at « = 1 (left column) and z = 10 (right column). 
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FIG. 8. Examples of two-component vortex solitons supported by the third-harmonic generation at s = —1: (a) a = 9, and 
(b) a = 13. 
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